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Abstract The behavior of an ultrasonic shot peening 
process is observed and analyzed by using a model of in- 
elastic hard spheres in a gravitational field that are flu- 
idized by a vibrating bottom wall (sonotrode) in a cylin- 
drical chamber. A marked heterogeneous distribution of 
impacts appears when the collision between the shot and 
the side wall becomes inelastic with constant dissipation. 
This effect is one order of magnitude larger than the sim- 
ple heterogeneity arising from boundary collision on the 
cylinder. Variable restitution coefficients bring the sim- 
ulation closer to the general observation and allows the 
investigation of peening regimes with changing shot den- 
sity. We compute within this model other physical quan- 
tities (impact velocities, impact angle, temperature and 
density profile) that are influenced by the number N of 
spheres. 



1 Introduction 

It is well known that the introduction of residual com- 
pressive stresses in metallic components leads to reduce 
fatigue strength Therefore, many engineering tech- 
niques involve surface treatment to allow either surface 
hardening (by e.g. nitriding or vapor deposition) or fa- 
tigue life improvement Q through laser shock peening 
or shot peening. For the latter, a high velocity stream of 
steel particles is projected at a material surface produc- 
ing at and below it compressive residual stresses with a 
peak value being reached at some depth below the sur- 
face ■ A particular mechanical treatment derived from 
conventional shot peening is called ultrasonic shot peen- 
ing, ultrasonic is a reference to the frequency of vibration 
of the sonotrode (see below) . It has received attention in 
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the recent years since it could be a promising tech- 
nique for obtaining surface treatment of metallic sur- 
faces. 

Here, a piezoelectric generator produces the vibra- 
tion of a sonotrode that projects upon contact steel shot 
in a chamber closed by a cover which is the sample to be 
peened. The shot is usually made of small steel particles 
whose diameter is between 1 and 3 mm and the fre- 
quency is about 20kHz. Several parameters can also be 
changed, allowing one to control of the overall shot veloc- 
ity and thus the shot peening intensity. Possible tuning 
parameters for optimizing the peening process include 
the shot diameter, the height of the chamber, the am- 
plitude and/or the frequency of the sonotrode. Basic ap- 
plications of this technology are found in automotive or 
aerospatial industry. 

If the performance of this process is closely related 
to the appropriate choice of parameters, it becomes nec- 
essary to understand how the peening intensity or the 
peening distribution on a given sample is affected by 
changes of mechanical or electrical characteristics of the 
system. Furthermore, as some of the physical quantities 
involved in the peening process are hard 1 to measure in a 
real-time experiment (velocity, acceleration,. . . ), mainly 
for safety reasons (the impact of steel beads are so strong 
that the chamber is a closed box), numerical simulation 
can be a powerful tool for investigating the influence of 
these quantities under various situations. 

In parallel, there has been great activity during the 
last ten years in the study of granular gases 0; B ■ ^ n 
particular, systems of vibro-fluidized glass beads in an 
cylinder has some similarities with the device used for 
the shot peening 0; [|| . The main difference with the ex- 
perimental setup used in ultrasonic shot peening is the 
absence of the cover and a lower frequency of vibration. 
In the latter series of experiments, it was shown that the 
inelastic sphere model provides an accurate description 
of microscopic quantities (local granular temper- 

ature, local mean velocity, local density, ...) which en- 

1 Here, it is a real-time experiment. 
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courages us to perform a molecular dynamics of inelastic 
hard spheres. 

The system is represented by a collection of inelas- 
tic hard spheres colliding with each other and with the 
boundaries (chamber, the shot, and sonotrode). For the 
sake of simplicity, collisions are first characterized by a 
constant normal coefficient of restitution. We perform an 
event-driven Molecular Dynamics that is as close as pos- 
sible to the experimental setup by using the geometrical 
features of the chamber, sonotrode and the cover. In or- 
der to obtain an improved description of the model, we 
also consider a model of inelastic hard spheres where the 
restitution coefficient depends on the relative velocity of 
the impact. 

Our results both theoretical and experimental show 
that the peening distribution on the sample is not homo- 
geneous. The heterogeneity of the peening distribution is 
strongly influenced by the value of the particle-side wall 
coefficient of restitution c w . This result goes far beyond 
the intuitive view that heterogeneity should simply re- 
sult from the boundary collisions on the side walls. The 
increased energy dissipation along the side walls favors 
particle accumulation thus increasing the gas (shot) den- 
sity on the border of the chamber. This leads to an in- 
crease of the impact frequency on the border of the sam- 
ple. Within the model, we compute impact velocity and 
impact angle of the shot and show also a changing behav- 
ior with the shot density, ranging from the dilute Knud- 
sen limit to a more dense situation where inter-particle 
collisions dominate. Both quantities display marked dif- 
ferences between the border and the center of the top 
wall (sample). 



2 Inelastic hard sphere model 

2.1 Simulation details 

We first consider the model close to the experimental 
setup (see below). The cylinder has a radius of R = 
35 mm which contains N = 200 hard spheres repre- 
senting the shot of diameter 3mm. The latter are sub- 
ject to a constant gravitational force. The energy is sup- 
plied by vibrating the bottom wall following a symmetric 
saw-tooth profile with amplitude A and period T which 
mimics in the simulation the sinusoidal profile of the 
sonotrode 2 . One should note that the choice of this pro- 
file has no major impact on the results 12], since the am- 
plitude of the harmonic n of the saw-tooth profile falls 
as odd n -2 . We mention also that even though the elec- 
trical excitation of the sonotrode is sinusoidal, because 

2 By using a saw-tooth profile the time of a collision be- 
tween a particle and the base is obtained analytically , 
whereas with a sinusoidal profile, the collision time is given 
by an implicit equation which requires a more expensive nu- 
merical computation. 



of the elastic deformation of the sonotrode, the velocity 
applied to the shot is certainly not purely sinusoidal 3 . 

The spheres collide inelastically and instantaneously 
with each other, with the cylindrical side walls, with the 
top wall and with the sonotrode. The corresponding con- 
stant coefficients of restitution are denoted c, c w , c& and 
Ct- The different collision rules are given by the following 
expressions: 

v i,r = v i,r - (1 + C t i))(Vi, r .r i>r )f i)T . (1) 
v 'i,z = v i,z - (1 + Cb){vi,z ~ vs) (2) 

v 'i,z = -C 1 + c t)vi, z (3) 
1 + c 

< ,3 = ± ~y~ I( v j ~ v *)- n l n ( 4 ) 

where the prime quantities denote the post-collisional 
quantities; v,* )7 . and f ^ is the unit position vector of par- 
ticle i are the velocity and the position of the particle i in 
the horizontal plane respectively and c w the normal coef- 
ficient of restitution for a collision between a sphere and 
the chamber. The particle-bottom (sonotrode) wall resti- 
tution coefficient Cb is first taken as unity which amounts 
to rescaling the amplitude of the vibration. v$ is the ver- 
tical velocity of the sonotrode ; v^ z is the vertical com- 
ponent of the velocity of particle i and ct the normal 
coefficient of restitution for a collision between a sphere 
and the cover; finally, a denote the velocities of i or 
j particle, n is the unit center-to-center vector between 
the colliding pair i and j and c is the normal coefficient 
of restitution for a sphere-sphere collision. One should 
note that there are two rules along the z axis (equations 

and ||3J|) depending on which wall the spheres collide: 
bottom (sonotrode) or top (sample). Between collisions, 
the spheres follow parabolic trajectories due to the con- 
stant gravitational field (viscous damping with the air 
contained in the chamber is neglected). 

The parameters of the model used in the simulation 
are obtained from experiment (Figure 0). The vibration 
frequency is 20kHz and the amplitude of the sonotrode 
is 25/jto. The chamber height is 40 mm. In the first part, 
we carry out most of the simulations with the restitution 
coefficients c = ct = 0.91 corresponding to the usual ex- 
perimental values for steel shot in the velocity range of 
interest [lflj and leaving as an adjustable parameter only 
the side wall restitution coefficient c w . This is in order to 
highlight the strong influence of the side-wall collisions 
on the impact heterogeneity. We stress however that the 
overall observed behavior of the model does not depend 
on the precise values given for ct and c and several addi- 
tional runs with different c t and c have been performed 
to check the robustness of our conclusions. The station- 
ary non-equilibrium state is achieved by a preliminary 
simulation of typically 5000 collisions per particle which 

3 There can be a difference in amplitude of the sonotrode 
on the border and on the center. 
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TOP WALL (sample) 

Chamber (H=40mm, R=35mm) 
Shot 



BOTTOM WALL (sonotrode; 
f=2()kHz, a=25 microns 



Fig. 1 Sketch of the experimental setup that is used as 
simulation box in the inelastic hard sphere model. 



corresponds to a peening time of about lms. The colli- 
sion time estimated as 1/u.s from Hertz theory EH is much 
lower. The statistical analysis of the quantities of inter- 
est has been accomplished for a total simulation time of 
26 s corresponding to 5 x 10 6 collisions. 



2.2 Velocity dependent restitution coefficient 

In a second part fsection ET^ . the present model is made 
more realistic by taking into account velocity dependent 
normal restitution coefficients that depend on the normal 
impact vel ocity . This dependence is rather well known 
El E E III El and was first reported in the 1920's. 
Plastification under high velocities (typically when v > 
5 m.s -1 ) has also been reported. In the high velocity 
limit, experimental measurements suggest |l4j a power- 
law behavior of the form: cn oc i> -1 / 4 whereas in the 
low velocity range, the deformations are supposed to be 
elastic and dissipation described by visco-elasticity [53 
l2lj . For the latter, it has been obtained j2£j a slightly 
different power-law which is like (1 — cjv) oc w" 1 / 5 . 

The behavior of restitution coefficients with respect 
to some easily measurable parameter is, however, a much 
deeper problem that can certainly not be encoded in the 
simple aforementioned power-laws. Experiments have in- 
deed shown that cn could depend on the sphere density 
p^ . the sphere diameter or the thickness of the impacted 
surface j2J| , or even the impact angle [3 • Beyond 
the details of a given material, studies on the restitu- 
tion coefficient all suggest a generic threshold between 
a regime at low velocity and low dissipation where cn 
depends weakly (2|| |2]| on the impact velocity, and a 
more dissipative regime induced by plasticity (or even 
fracturing) pl : ! ^ . 

Recently, simulations of inelastic hard spheres have 
been performed using variable restitution coefficients [23| 
and have shown the necessity of the latter to accurately 
describe experiments. Specifically, pressure effects as a 
function of the density of spheres could be recovered 
by simulation for a dilute and dense vibrated granular 



Impacted material 


i 


CO 


Wo [cm/s] 


Sonotrode (titanium) 


b 


0.91 


1.2 


Spheres (steel) 


s 


0.91 


1.2 


Sample (aluminum) 


t 


0.6 


0.12 


Side walls (aluminum) 


w 


0.6 


0.12 



Table 1 Parameters for the inelastic hard sphere model with 
variable restitution coefficient. 



medium. It has been shown also that the unphysical clus- 
tering tendency was reduced with the use of velocity de- 
pendent restitution coefficients. 

In section IIIB, we will use a threshold model for the 
normal restitution coefficient defined by: 




v < v. 



0' 



V > Vq 



(5) 



where v q is a threshold velocity, Cq is the constant normal 
restitution coefficient at low velocity and i = b,t,w,s fol- 
lowing the nature of the impacted surface (bottom, top, 
wall, spheres). Parameters are given in Table I. One ex- 
pects indeed that softer materials such as the aluminum 
side walls or the sample will have a lower threshold veloc- 
ity and a lower Co as compared to the titanium sonotrode 
or the steel spheres. 

Finally, we also use the simplest possible model to 
take into account the transverse dissipation that leads 
to a tangential restitution coefficient ct- To our knowl- 
edge, only a very few studies have been reported on the 
subject (see however |^|12)- The total loss of trans- 
lational kinetic energy can usually be described by the 
total restitution coefficient c = [c N cos 2 9 + c\ sin 2 0] 1 / 2 
where 9 is the impact angle. Conservation of impulse 
and momentum and an additional condition of rolling 
prior to departure from the impacted surface leads to a 
value of ct = 5/7 (the factor 5/7 comes from the mo- 
mentum of intertia). Here, it is assumed that the loss 
in kinetic energy mostly arises from the sphere rotation 
during the impact. One may also assume that slip con- 
tinues throughout contact which will in this case lead 
to: ct = 1 — /i(l + cn) cot 9. But this would define an 
additional parameter fi corresponding to the Coulom- 
bic friction coefficient of the impacted surface. Experi- 
mental measurements of steel spheres bouncing on flat 
aluminum plates show [24j that the constant value of 
ct = 5/7 is mostly valid at small impact angles and un- 
der certain conditions up to 9 ~ 55°. However, a more 
detailed analysis that should include the deformations 
and velocities associated with the elastic deformations 
of the surfaces is clearly beyond the scope and objec- 
tives of this paper. Since we are handling instantaneous 
impacts, effects of friction or peculiar material properties 
can only be taken into account via an effective velocity 
dependence of the restitution coefficients. Furthermore, 
we stress that this will not affect the general observed 
behavior with shot density. 
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3 Results 

3.1 Origin of the impact heterogeneity 

Figure [21 shows the different impact profiles that appear 
on the top wall of the chamber after 1 s simulation time 
for two different values of the side wall restitution coef- 
ficient c w . For c w — c — ct — 0.91, one has an almost 
homogeneous distribution of impacts (Fig whereas 
heterogeneity sets in when c w is lowered to 0.20 (Figure 
The present results have to be contrasted with the 
observed profiles obtained on the impacted aluminium 
sample after 1 s or on the sonotrode. It clearly shows 
that the inelastic sphere collision with the side wall are 
relevant for understanding the heterogeneous shot peen- 
ing which is manifested by an increased impact number 
on the border of the sample. Nonelastic collisions on the 
side wall originate the impact profile that is experimen- 
tally observed (e.g. on the sonotrode (bottom wall), Fig 
Ht) and recovered from the simulation. The same pro- 
files are obtained for the top wall with a similar impact 
frequency (see below). 

The influence of the sphere- wall coefficient of restitu- 
tion is quantitatively observed by monitoring the impact 
frequency per surface unit Af with respect to the radius R 
of the chamber (Fig.|3J). It results that when the shot has 
a pure elastic collision on the side walls {c w = 1.0), then 
the number of surface impacts slightly increases from the 
center to the radius of the chamber, whose values are in 
the range Af = 0.6 — 0.8 impact. mm~ 2 .s^ 1 . When c w 
decreases, a similar behavior is observed when the radial 
distance is lower than a typical value i?o — 30 mm. For 
R > Ro significant differences appear, e.g., for c w — 0.60, 
the number of surface impacts is multiplied by almost a 
factor 10 between the center and the border of the cover 
and this increase is even more dramatic for the ultimate 
value of c w = 0.20. 

Intuitively, one expects that the impact distribution 
should be heterogeneous even in the case of elastic side 
wall collisions because of the lateral bounce on the cylin- 
der wall. The simulation shows that this effect clearly 
seen from Fig. |3| when c w — 1.0 is weak as compared to 
the effect arising from the decrease of c w . Boundary col- 
lisions produce an increase of the surface impacts from 
~ 0.5 mm~ 2 .s^ 1 in the center to about 0.8 mm~ 2 .s 
on the border for elastic side wall collisions whereas the 
effect of lower side wall restitution coefficients (e.g. c w = 
0.6) leads to Af — 1.8 mm~ 2 .s~ l on the border of the top 
wall, compared to VV = 0.25 mm~ 2 .s _1 on the center. 
The present results do not depend crucially on the val- 
ues taken for the other restitution coefficients as similar 
trends for the surface impact frequency are obtained for 
lowered c and Cj. For instance, when c = ct — c w — 0.6, 
the trend with R observed is rather close to the one dis- 
played in Fig. |3]with c w = 0.6, except that Af ranges 
now from 0.1 to 1.7 rnmT 2 .s~ x between the center and 
the border of the sample. 




Fig. 2 Distribution of impacts on the top wall (sample), a) 
c w — 0.91. b) c w = 0.20. c) Sonotrode after several hours of 
use. See text for details. 
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Radius [mm] 

Fig. 3 Surface impact frequency with respect to the radial 
distance R for different side wall restitution coefficients c w 
(symbols). The broken line corresponds to the pure elastic 
side wall collision (c w = f). The distance Ro is used below 
(see text). 

Beyond the numerical details, the origin of the het- 
erogeneity becomes clear. With increased dissipation on 
the side walls, the spheres have a reduced velocity and 
are "adsorbed" on the side walls with an upward heli- 
coidal trajectory arising from the impulse of the sonotrode 
As a result, the density and the granular temperature 
(kinetic energy per particle) of the granular gas appear 
to be also strongly influenced by the dissipation. 

3.2 Comparison with observation 

Aluminum samples have been peened during 1 sec and 
the impact distribution displays a weak heterogeneity. 
The average roughness between the border and the cen- 
ter of the sample is respectively 4 fj,m and 7.4 \xm. In gen- 
eral, the obtained roughness reveals the degree of impact 
on a sample. It provides therefore an indirect evidence 
that more shot has been impacting the border. The het- 
erogeneity can be quantified in a manner similar to FigOU 
by sampling circularly the frequency of impacts Af per 
unit surface with respect to the radius of the sample. Af 
varies from 0.65 at the center of the sample to 0.95 on the 
border. The same tendency is obtained with rectangular 
sampling, i.e. one obtains Af = 1.25 ± 0.12 at the center 
of the sample, and Af = 1.66 ± 0.39 at the border. The 
presence of polymer adhesive stripes on the side walls 
that induce an increased dissipation, leads to a lowering 
of the respective impact frequencies, in harmony with 
the observed trends displayed in Fig. [3] Refinement of 
the measurements is currently under consideration. 

3.3 Selective shot peening 

Once the origin of the heterogeneity is identified, we in- 
vestigate the effect of the sphere density (or sphere num- 




Impact angle [deg] 



Fig. 4 Impact angle distribution on the top wall for different 
number of spheres N. All distributions are normalized to one. 



ber Af) in the chamber on the peening statistics. Within 
the inelastic hard sphere model using now the " realis- 
tic" restitution coefficients of Eq. (JSJ) , one has indeed the 
luxury to investigate features appearing with the pro- 
gressive jamming of the system, starting from the very 
dilute limit, and characterize the deeper origin of the 
obtained profiles. 

Figure 01 shows the impact angle distribution on the 
top wall (the sample) for various number of spheres. 
One can observe that the nature of the impacts can be 
very different following the very dilute (N = 50, den- 
sity r]o = 3.25 x 10~ 4 mm -3 ) or more dense (N = 1000, 
r/o = 65.0 x 10~ 4 mm -3 ) situation. At low density, the 
impact is almost normal with a very sharp distribution 
centered around the impact angle 8 = 0°. Here one sees 
that the spheres will mostly bounce back and forth be- 
tween the top and the bottom walls with a rather small 
number of inter-particle collisions. As a result, the sur- 
face impact frequency Af with respect to the radius (Fig. 

is rather flat and starts only to grow close to the side 
walls. With increasing N, these simple (mostly linear) 
trajectories tend to disappear as more and more sphere- 
sphere collisions are now involved. Finally, the distribu- 
tion becomes very broad at high N and centered around 
9 = 35° and there is not much difference between the 
system with N = 600 and N = 1000 spheres. Addi- 
tional spheres do not change the obtained distribution. 
This shows also that in the more dense situation, normal 
impacts are very rare on the top wall, i.e. the probabil- 
ity of having 9 — is almost zero. Noteworthy is the 
system with N = 200 spheres which displays an angu- 
lar distribution that contains a reminiscent signature of 
the very dilute limit, i.e. showing a shallow peak around 
9 = 5° which disappears when N is increased from 200 to 
400. It suggests that in this intermediate situation, some 
spheres succeed in moving upwards through the granular 
gas without any sphere-sphere collision. 
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Fig. 5 Surface impact frequency J\f on the top wall as a 
function of the radius of the chamber for different number of 
spheres. 




Velocity [m/s] 



Fig. 6 Vertical impact velocity distribution on the top wall 
for different number of spheres N. All distributions are nor- 
malized to one. 



Figurc[S]sliows an additional interesting feature which 
is the occurrence of a local order produced by the accu- 
mulation of the spheres on the side walls. With increas- 
ing N, more and more spheres are trapped on the side 
walls, similarly to the result of the constant restitution 
coefficient model (Figure |3J). This produces an increased 
jamming in the vicinity of the side wall which does not 
allow other arriving spheres to reach it. As a result, these 
spheres will stay at a distance of the order of the diame- 
ter a of the spheres. This is reflected in the quantity M 
by an impact frequency peak at about R = 30 mm and 
even a secondary peak for high densities (N — 600 and 
N = 1000) between the latter value and the border of 
the sample. 

Finally, we note that for a large number of spheres 
(N = 1000), the surface impact frequency decreases from 
the center of the sample to the first impact frequency 
peak at R — 30 mm. With the increased number of 
spheres, the impact velocity distribution (Fig. refines 



10 20 30 

Radius [mm] 

Fig. 7 Surface impact frequency M on the top (solid lines, 
same as Fig. [3J and the bottom (broken lines) walls as a 
function of the radius of the chamber. 



and converges to a Maxwell-Boltzmann-like distribution 
that can be fitted by v z exp(— u|) with a mean velocity 
that is of about 0.8 m/s and which is close to the mea- 
sured velocities for this kind of system Note that 
for a small number of spheres, the impact velocity dis- 
tribution is broad and ranges from 1 m/s up to large 
velocities of about 6 m/s corresponding to spheres that 
have been optimally accelerated by the sonotrode. Once 
the system densities, the dissipation due to sphere-sphere 
collisions lowers the overall impact velocities. 

A comparison between the surface impact frequency 
of the top and the bottom walls (Figure [7J| shows that 
the impact regime can be rather different but still con- 
sistent with previous findings. For a very small number 
of spheres (N — 50), the result on the surface impact 
frequency is consistent with our previous findings, i.e. Af 
is identical between the top and the bottom, in agree- 
ment with the conclusion drawn from the impact an- 
gle distribution (Fig. QJ, i.e. quasi- normal trajectories 
for the spheres which collide almost only between bot- 
tom and top. For an increased number of spheres, dif- 
ferences emerge which can become very pronounced as 
suggested by the value of the impact frequency on the 
center of sample for N = 1000. For the top wall, it is 
found M = 0.50 mm 2 .s _1 whereas M — 0.26 mm 2 .s _1 
for the bottom wall, i.e. there is almost a factor of two 
between the two colliding walls. In order to infer the ori- 
gin of the peening difference of the two surfaces, we have 
computed the mean velocity field for the system of in- 
elastic spheres in the (i?, z) plane. The field is averaged 
over the azimuthal angle and 5 x 10 6 collisions. Figure|H| 
shows the velocity field for the corresponding number of 
spheres used in Fig. [7| For N = 200, a toroidal convec- 
tion roll is clearly present in which the particles flow, on 
average, up from the border and down the center. This 
kind of convection roll has been found both in simula- 
tions and experiment for open vibrated granular media 
§ El at vibration frequencies of 50 Hz. The convec- 
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Fig. 8 Mean velocity field of the inelastic hard spheres in the 
(r,z) plane for three different numbers of spheres, a) N=50, 
b) N=200, c) N=1000. 



40 



30 



S 20 



10 




- oo 

II 



o%Vo 



10 20 

Radius [ram] 

Fig. 9 A snapshot of the N = 1000 shot inside the chamber 
showing the density difference close to the top and close to 
the bottom wall. An average has been performed over the 
azimutal angle, which explains observed overlaps discs. For 
clarity, the size of the spheres has been reduced. The broken 
lines serve to de fine the regions used in the discussion (see 
text) and in Fig.lTTil 



tion roll is maintained (but weaker by about an order 
of magnitude in intensity) for a lower number (N = 50) 
of spheres. However, one sees that for N = 1000, the 
convection roll breaks up and spheres flow from the side 
wall either upwards to hit the top of the chamber, or 
downwards to the sonotrode. 

This tends to separate the chamber into two parts 
(Fig. A first part (at z > z$ with zq ~ 20 mm) 
that connects to the top wall where the density is large 
(rj = 8.12 x 10~ 3 mm~ 3 and a packing fraction of 0.115). 
On the other hand, the lower part of the chamber cor- 
responds to a much more dilute situation (j) = 4.87 x 
10~ 3 mm~ 3 and a packing fraction of 0.069). Conse- 
quently, the upper impact angle distribution is radically 
different. The impact angle of the spheres bouncing on 
the more dense media at z > zq is very close to 9 = 0°. 
This contrasts with the impact angle distribution of the 
upper part (Fig. fTTTf) . 



4 Summary and conclusion 

Simulations on inelastic hard sphere models using con- 
stant restitution coefficients clearly show that the inelas- 
ticity of the side-wall collisions plays a key role in the 
impact profile of the bottom and the top walls of a peen- 
ing chamber. With increased dissipation, an increased 
heterogeneity of the impacts arising from the accumula- 
tion of the spheres on the side walls is found. The cor- 
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Fig. 10 Upper impact angle distribution for spheres be- 
longing to region I (z > Zo, solid line, same as Fig. |1J and 
region II (z < z , R < R broken line) for N = 1000. Both 
distributions are normalized to one. 



responding surface frequency of impacts shows that the 
latter effect is one order of magnitude larger than the 
heterogeneity produced simple by oblique collisions aris- 
ing from the side walls. 

A model using variable restitution coefficients per- 
mits us to study in more detail the effect of the shot 
density (or number N of spheres) in the chamber. It 
shows that different peening regimes on the top wall take 
place with changing N that range from normal impacts 
for dilute granular gases, to oblique impacts with a well- 
defined mean impact velocity. Densfication close to the 
side walls produces the occurrence of a local order at 
a distance of about the sphere diameter from the side 
walls. 

These results suggest that elastic control of the side 
wall and a careful selection of the shot density will per- 
mit us to tune peening regimes for the shot and allow 
various kinds of surface treatments. The high density 
observed in the vicinity of the side walls is associated 
with a downwards helicoidal trajectory of the spheres. 
We believe that the use of a very rough surface on the 
side wall could lead to the reinjection of the spheres in 
the bulk. Further consideration in this direction along 
with deeper experimental characterization of the hetero- 
geneity is currently under consideration. 
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